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ABSTRACT 

We present a combined analysis of mass estimates in the central cores of galaxy clus- 
ters from the strong lensing, the X-ray measurements and the universal density profile 
(NFW). Special attention is paid to the questions (l)whether the previously claimed 
mass discrepancy between the strong lensing and X-ray measurements is associated 
with the presence of cooling/non-cooling flows, (2)whether the cusped NFW den- 
sity model can provide a consistent cluster mass with the strong lensing result and 
^vq | (3)whether a non-zero cosmological constant can be of any help to reducing the strong 

lensing - X-ray mass ratios. We analyse a sample of 26 arc-like images among 21 clus- 
\£) • ters,the X-ray data of which are available in archive. The X-ray and NFW cluster 

\ masses are obtained by assuming that the intracluster gas is isothermal and in hydro- 

• static equilibrium with the underlying gravitational potential of the clusters. 

A statistical comparison of these three mass estimates reveals that the mass dis- 
crepancies for all the events are well within a factor of 2, if X-ray measurement uncer- 
fS tainties are included. In particular, we confirm the result of Allen that the larger mass 

discrepancy is only detected in the intermediate cooling, especially non-cooling flow 
clusters, thus attributing the mass discrepancy to the local dynamical activities in the 
central regions. We show that the NFW profile yields a consistent cluster mass with 
the conventional X-ray measurement, which is interpreted as the consequence of the 
common working hypothesis behind the two methods. Any difference between these 
two models must occur at even smaller radii (e.g. within the arc-like images) or at 
large radii. It appears that the introduction of the cusped density profile as the dark 
matter distribution of clusters cannot raise the cluster masses enclosed within arc-like 
' images. Finally, a non-zero cosmological constant is able to moderately reduce the 

mass ratios of m\ oria to m xray . These results, together with the excellent agreement be- 
tween the X-ray, optical and weak lensing determined cluster masses on scales greater 
than the X-ray core sizes found in the early work, indicate that the mass discrepancy 
between strong lensing and other methods in the intermediate cooling and non-cooling 
flow clusters is likely to have arisen from both the oversimplification of lensing model 
and the inappropriate application of isothermality and equilibrium hypothesis in the 
central regions of clusters where the local dynamical activities make a non-negligible 
contribution to both mass estimates. 

Key words: clusters: general - dark matter - galaxies: gravitational lensing - X-rays: 
galaxies. 
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1 INTRODUCTION 

An accurate estimate of the total gravitating masses of clus- 
ters of galaxies is crucial for determinations of the mass- 
to- light ratios M/L and the baryon fractions /b of clus- 
ters, while the latter play a potentially important role in 
the 'direct' measurement of the mean mass density of the 



Universe, Qm- Four different techniques have been available 
today to determine the masses of clusters: the optical mea- 
surements of distribution and velocity dispersion of cluster 
galaxies, the X-ray measurements of intracluster gas and 
its temperature, the gravitationally distorted images of dis- 
tant galaxies behind clusters, and the numerical simulations 
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of formation and evolution of clusters. Because traditional 
cluster mass estimators using optical/X-ray observations of 
galaxies/intracluster gas rely strongly upon the assumption 
of hydrostatic equilibrium, the gravitational lensing serves 
as a powerful and efficient tool for both estimating cluster 
masses and testing the accuracy and reliability of traditional 
methods. It turns out that there is good agreement between 
the gravitational lensing, X-ray and optical determined clus- 
ter masses on scales larger than the X-ray core radii, within 
which the X-ray method is likely to underestimate cluster 
masses by a factor of 2-4 ( Wu 1994; Miralda-Escude & Babul 
1995; Wu & Fang 1997; Allen 1998; Wu et al. 1998 and ref- 
erences therein). While a number of mechanisms have been 
suggested for the reported mass discrepancy (Loeb & Mao 
1994; Miralda-Escude & Babul 1995; Allen 1998), a satisfac- 
tory explanation has not yet been achieved. It is generally 
believed that such a mass discrepancy has probably arisen 
from either the oversimplification of strong lensing model 
for the central mass distributions of clusters (Bartelmann 
& Steinmetz 1996) or the inappropriate application of hy- 
drostatic equilibrium hypothesis in the central regions of 
clusters (Wu 1994; Wu & Fang 1997). 

In this paper we intend to have a close examination 
of the issue as to whether the reported mass discrepancy 
is associated with the following three factors: (l)the pres- 
ence of cooling/non-cooling flows in clusters, (2)the cusped 
dark matter profile, i.e. the so-called universal density pro- 
file (Navarro, Frenk & White 1995; hereafter NFW), as the 
mass distribution of clusters, and (3)a non-zero cosmological 
constant (£Ia)- 

Our reconsideration of the influence of cooling/non- 
cooling flows on the X-ray mass estimate of clusters is mo- 
tivated by the pioneering work of Allen (1998), who has 
demonstrated that there is an excellent agreement between 
the strong lensing and X-ray determined cluster masses for 
cooling-flow clusters while the mass discrepancy is only de- 
tected among the non-cooling flow ones owning to the signif- 
icant offsets between the X-ray and lensing centers. If this is 
the case, the presence or absence of cooling-flows can be used 
as an indicator of whether or not clusters have reached the 
state of dynamical relaxation. Meanwhile, the mass discrep- 
ancy between the strong lensing and X-ray measurements 
can be attributed to the local dynamical activities of clus- 
ters where the hydrostatic equilibrium assumption becomes 
invalid. However, such a scenario has been argued recently 
by Lewis et al. (1999) based on an analysis of 14 CNOC 
clusters. They have shown that the systematic effects own- 
ing to cooling-flows, non-equilibrium systems and tempera- 
ture gradients on the average ratio of the 'dynamical' cluster 
masses to the X-ray masses do not exceed 15-20 per cent. 
Therefore, further work is needed to clarify the issue. 

The original purpose of this investigation was to study 
whether the cusped NFW profile as the dark matter distri- 
bution of clusters can resolve the discrepancy between the 
strong lensing and X-ray mass estimates. First, the central 
singularity in the NFW profile is consistent with the early 
claim, based on the modeling and statistics of strongly dis- 
torted images of background galaxies by massive clusters 
(Hammer 1991; Wu & Hammer 1993; Grossman & Saha 
1994), that dark matter profiles are sharply peaked towards 
cluster centers. Secondly, except for the small core radius 
there is a striking similarity between the distribution of in- 



tracluster gas tracing the dark halo of the NFW potential 
and the conventional /3 model (Makino, Sasaki & Suto 1998). 
Thirdly, many high-resolution simulations of structure for- 
mation have shown that the NFW profile is independent of 
mass, initial density fluctuation or cosmology (e.g. Cole & 
Lacey 1996; Navarro, Frenk & White 1997; Eke, Navarro 
& Frenk 1998). In particular, Makino & Asano (1999) have 
recently applied the NFW profile to three lensing clusters, 
A2163, A2218 and RXJ1347, and compared the NFW profile 
and strong lensing derived cluster masses. Their study indi- 
cates that the steeper cusps in the NFW profile can indeed 
reduce the mass discrepancy between the X-ray and gravita- 
tional lensing measurements. On the other hand, the NFW 
profile can also recover the observed surface number density 
of cluster galaxies (e.g. Carlberg, Yee & Ellingson 1997), 
from which one can fix the two free parameters in the NFW 
profile and then work out the total dynamical masses of the 
clusters. Using the 14 CNOC distant clusters, Lewis et al. 
(1999) found no systematic bias between the (NFW) dynam- 
ical and X-ray methods. Consequently, the dynamical mass 
derived from the cusped NFW profile cannot be reconciled 
with the gravitating mass from strong lensing. This implies 
that either the NFW profile may not be extrapolated to 
smaller radii (~ 10 kpc) or the equilibrium hypothesis may 
break down in the central regions of clusters, provided that 
the current strong lensing model gives a reliable estimate 
of the cluster masses enclosed with the arc radii. Alterna- 
tively, the giant arc statistics, based on the NFW profile as 
the underlying gravitational potential and the X-ray mea- 
surements of the gas distribution, can marginally reproduce 
the observed number of giant arcs but requires unreasonably 
high X-ray temperatures for some clusters (Molikawa et al. 
1999). Thus, the detailed modeling of strong lensing includ- 
ing the 'irregularity' in the mass distribution of a cluster is 
once again advocated. Taking these results as a whole, we 
feel that at least one of the following working hypotheses 
should be abandoned: (1) a simple estimate of cluster mass 
from strong lensing, (2) the hydrostatic equilibrium in the 
central region of a cluster, and (3) the NFW profile as the 
total mass distribution. In order to stress the point, more 
examples will be provided in this paper and a systematic 
comparison among the cluster masses given by the strong 
lensing, the NFW profile and the X-ray measurements will 
be made. 



Another question we would like to address is how our 
estimates of the strong lensing and X-ray cluster masses are 
affected by a non-zero cosmological constant Qa- This arises 
because the lensing mass miens is proportional to DdD s / Dds 
while the X-ray mass m xr ay goes as Dd, where Dd, D a and 
Dds are the angular diameter distances to the cluster, to 
the background galaxy, and from the cluster to the galaxy, 
respectively. As a result, the mass ratio mi en s/wixray de- 
pends on the cosmic density parameter !1m and J1a through 
D B /Dds- In particular, mi cns /m xray becomes smaller in an 
Qa dominated universe than in an Qm dominated one. This 
opens a new possibility to reduce the large discrepancy be- 
tween miens and m X ray We will compare the lensing and 
X-ray cluster masses in a flat cosmological model with and 
without the cosmological constant. Throughout this paper 
we assume Ho = 50 km s _1 Mpc -1 . 
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2 CLUSTER MASS ESTIMATES 
2.1 Strong gravitational lensing 

Gravitational lensing furnishes a simple yet efficient way to 
measure the projected cluster mass along the line of sight, 
ft is believed that a simple spherical lensing model provides 
rather a good estimate of the projected cluster mass within 
the position (r arc ) of arc-like image, which reads 



TOlc 



(< 

Tare J — TTTarc^-'c 



(1) 



where S cr it = {c 2 / AnG)(D s / D^Dds) is the critical surface 
mass density. The above equation is actually the lensing 
equation for a cluster lens of spherical mass distribution 
with a negligibly small alignment parameter for the distant 
galaxy with respect to r arc . We adopt the value mi cns (< r arc ) 
from the detailed modeling of arclike images where available 
from literature. 



2.2 X-ray measurements 

Assuming that the diffuse X-ray emitting gas in a cluster is 
isothermal and in hydrostatic equilibrium with the underly- 
ing gravitational potential of the cluster, we have 

GM(r) _ kT dlnn gas (r) 



dr 



(2) 



where M(r) is the total cluster mass enclosed within ra- 
dius r, T and n gas (r) are the gas temperature and num- 
ber density, respectively, and fi = 0.585 denotes the mean 
molecular weight. If the spatial distribution of intracluster 
gas is described by the conventional f3 model, n gas (r) = 
n gas (0)(l + r 2 /r 2 )" 3 ' 3 '' 2 where r c is the core radius of the 
X-ray gas profile, we can easily get the dynamical mass dis- 
tribution from equation (2). In order to compare with the 
gravitational lensing result, we use the projected X-ray clus- 
ter mass within radius r (Wu 1994) 



m xray = 1.13 x 10 13 /3 flt m(r) 



kT 



0.1 Mpc J VI keV 



M , (3) 



where 
m(r) = 



(R/r c ) 2 + l 



and R is the physical radius of the cluster and will be taken 
to be R — 3 Mpc in the actual computation. Our conclusion 
is unaffected by this choice. 

2.3 The universal density profile 

The virialized dark matter halo follows (NFW) 

(4) 



(r/r.)(l + r/r.)"' 

where p B and r B are the characteristic density and length, 
respectively. The projected cluster mass within radius r from 
equation (4) is simply 

„3 



In + 

l n .J- j r " 



4-wp s r; 

_ m rs+\/ rf-r 2 



r < r s ; 



arctan ^ T r * , r > r s . 



(5) 



In principle, the two parameters p s and r s can be determined 
from the observed X-ray surface brightness of intracluster 
gas or the surface number density of galaxies, in combina- 
tion with the measurement of X-ray temperature or velocity 
dispersion of galaxies. We will work with the X-ray data in 
this paper. 

Inserting the total mass M(r) deduced from the NFW 
profile into equation (2) will result in an analytic form of 
gas profile (Makino et al. 1998): 



s (x) = n gas (0)e "(1 + x 



a/x 



(6) 



in which x = r/r s and a = ATvGnm v p s rl /kT '. Since 
(0)e a , we introduce a background sub- 
tracted gas number density h(x) — n gas (x) — n gas (oo), which 
reads 



ft(«) = ^l[(l + «) a/ "-l] 



(7) 



This avoids the divergence in the computation of X-ray sur- 
face brightness below. Recall that an arbitrary cut-off radius 
was used by Makino & Asano (1999) for the same purpose. 
The X-ray surface brightness in the scenario of the optically 
thin and isothermal plasma emission is 



S x (9) oc 



h 2 (x)xdx 



'D d e/r s y/x 2 - D A 9/r s 
A straightforward computation yields 



(8) 



<* fe°Je s 



\ Of. fx 



ln(l+x)l 



dx, 



(9) 



in which ^ s = r s /Dd- Our task is thus reduced to finding out 
the two parameters a and r s from the observed X-ray surface 
brightness S x and gas temperatures T of clusters, and then 
to obtain the projected cluster masses from equation (5) . 



3 CLUSTER SAMPLE AND RESULTS 

21 clusters (Table 1) are selected from the Strong Lensing 
Cluster Sample of Wu et al. (1998) by requiring that both 
the X-ray temperature and the surface brightness profiles are 
well measured with ROSAT and/or ASCA. We have not in- 
cluded those clusters whose temperatures are estimated by 
indirect methods such as the X-ray luminosity-temperature 
correlation and the velocity dispersion-temperature correla- 
tion. These 21 clusters have a mean redshift of (zd) = 0.30, 
and contain 26 arc- like images. For the 11 arcs that have 
no redshift information, we estimate their lensing masses 

^■lens 

by assuming the mean redshifts of (z s ) = 0.8 and 2.0, 
respectively. We calculate the values of mi cns (< r arc ) for 
two cosmological models: (1) J7m = 1 and Qa = and (2) 
£Im = 0.3 and = 0.7. The resultant cluster masses en- 
closed within the positions of 26 arc-like images are listed in 
Table 1. We have not provided the uncertainties in rrii ens no 
matter whether they are derived from the detailed modeling 
of arcs or from the simple spherical assumption. In princi- 
ple, there will be no uncertainties associated with mi cns in 
the simple spherical lensing model if redshifts of the lensing 
clusters and of the arc-like images are reliably determined. 
However, the presence of substructures and asymmetrical 
mass distributions in clusters may complicate this simple 



© 0000 RAS, MNRAS 000, 000-000 



4 Wu 



cluster mass estimate, which will be the major source of un- 
certainties in miens. The absence of the secondary arc- like 
images in most of the arc-cluster systems should be a strong 
argument against the spherical mass distribution in the cen- 
tral regions of clusters. 

The X-ray surface brightness profiles S x for our 21 clus- 
ters are available in archive (ROSAT and ASCA), which 
have been analysed by a number of authors for different 
purposes. Because the accuracy of our determination of the 
dynamical cluster mass from X-ray measurement is closely 
connected to the issue as to how precisely we can describe 
the observed X-ray surface brightness S x , two conventional 
methods are employed in the fitting of S x (e.g. Neumann 
& Arnaud 1999): (1) a ft model fit to the entire X-ray sur- 
face brightness; (2) a j3 model fit by excising the central 
region until an improvement of \ 2 can be achieved. The 
latter is usually applied to the cooling-flow cluster where a 
sharp peak is seen in the X-ray emission concentrated in the 
core of the cluster. In Table 2 we list the best-fitting pa- 
rameters (3 and r c from the first fitting method, while the 
results obtained from the second method are presented in 
Table 3 for 14 clusters where available in literature. Here 
we make no attempt as far as possible to extrapolate the 
original work. For all the clusters we take the temperature 
data from the literature (see Wu, Xue & Fang 1999; and 
references therein). In order to explicitly demonstrate the 
possible effect of cooling/non-cooling flows on the mass ra- 
tios of miens /m xray , the sample is classified into the massive 
cooling-flow (MC) clusters, the intermediate cooling-flow 
(IC) clusters and non-cooling flow (NC) clusters in terms 
of their cooling times (see Allen 1998). Roughly speaking, 
the MC, IC and NC clusters in our sample have average 
mass deposition rates of <; lOOOM© yr _1 , ~ lOOM© yr _1 
and ^ lOM© yr~ , respectively. The projected dynamical 
masses (m xray ) of clusters within the arc positions are com- 
puted according to equation (3). Finally, the mass ratio of 
miens to m xray is obtained for each case and listed in Ta- 
bles 2 and 3, in which the error bars reflect the combined 
measurement uncertainties of T, r c and f3. 



In principle, we can get the two parameters a (or p B ) 
and r s in the NFW profile by fitting the observed S x to the 
predicted form of equation (9) for all the cases. However, 
the actual operation turns out to be difficult for some clus- 
ters owning to the instrumental PSF of the X-ray telescopes 
(ASCA, ROSAT). So, for eight clusters in our sample we 
estimate the parameters a and r s using the empirical for- 
mula (Ettori & Fabian 1999), a = 14.34/3 and r s = 3.17r c . 
For the three clusters (A2219, MS1358 and MS2137), we 
take the best-fitting values from Ettori & Fabian (1999) who 
have excluded the cooling-flow regions in their fitting. For 
the rest clusters, we adopt the results of Wu & Xue (2000) , 
which are obtained by the \ 2 fit of the entire data points 
of S x from the Mohr, Mathiesen & Evard (1999) sample to 
the theoretical predictions (equation 9). The parameters, a 
and r s , for all the clusters of our sample are listed in Table 
4, together with the derived cluster masses within the po- 
sitions of arcs and their ratios to the corresponding masses 



from strong lensing for a cosmological model of Qm = 1 and 
n A = 0. 

We first compare in Fig.l the projected cluster masses 
within arc-like images determined from strong lensing and 
X-ray measurements, using the data in Table 2, where the 
values of mie ns for z B = 2 are chosen if redshifts of the arcs 
remain unknown. It appears that the ratios of miens to m xray 
display rather large dispersions for many clusters, which is 
mainly due to the fact that the X-ray surface brightness 
profiles of these clusters (e.g. A370, CL0500, etc.) have not 
been well constrained. Yet, the mass ratios of all the data 
points essentially satisfy micns/m xray < 2 when their un- 
certainties are taken into account. Namely, the lensing and 
X-ray determined cluster masses are still consistent with 
each other within a factor of 2, and there is no large dis- 
crepancy between these two mass estimates. Additionally, 
all the clusters which exhibit a relatively large mass ratio 
of miens /m xray > 2, if the error bars are neglected for the 
moment, are the NC or IC clusters, in contrast with the MS 
clusters whose mass ratios mi cns /m xray are roughly consis- 
tent with unity. We have thus confirmed the results of Allen 
(1998). Alternatively, the introduction of a non-zero cosmo- 
logical constant does help to reduce, though moderately, the 
mass discrepancy between miens and m xray . 

Next we demonstrate how our X-ray mass measure- 
ments within r arc are affected by the two different fittings of 
the X-ray surface brightness profiles: whether or not the ex- 
cess X-ray emission in the central cores is excluded. We dis- 
play in Fig.2 the projected X-ray masses within 17 arc-like 
images obtained from these two methods for a cosmological 
model of Qm = 1 and J7a = 0. It is immediate that these two 
methods provide essentially a consistent cluster mass at the 
central core. Consequently, the mass ratios mi CIls /m xray will 
remain approximately the same as those in Table 2. even if 
the goodness of a (5 model fit is improved by omitting some 
data points of S x in the central regions. 

Finally, we compare the projected cluster masses within 
arc-like images derived from the strong lensing, the X-ray 
measurement based ona/3 model and the NFW profile. Fig. 3 
shows the results of itinfw(< J" arc ) versus m xray (< r arc ) for 
our sample of 26 arcs among 21 clusters. Regardless of the 
large error bars, the cluster masses given by the NFW profile 
are in good agreement with the values from the X-ray anal- 
ysis. Indeed, this result contrasts with our initial yet naive 
speculation that the cusped NFW profile at the cluster cen- 
ter can give rise to a very different cluster mass enclosed 
within the arc-like image. On the other hand, such a situ- 
ation is not surprising because both the X-ray and NFW 
cluster mass estimates are based on the measurements of X- 
ray surface brightness profiles of clusters and the assumption 
that the X-ray emitting gas is isothermal and in hydrostatic 
equilibrium with the underlying gravitational potential of 
the clusters. Recall that the NFW profile can almost recover 
the conventional j3 model for the observed X-ray surface 
brightness profiles of clusters. This explains the good agree- 
ment between these two mass estimates. The similar result 
has been reported by Lewis et al (1999) if the surface num- 
ber density of galaxies and their velocity dispersion, instead 
of the X-ray surface brightness profile of intracluster gas and 
its temperature, are used to calibrate the NFW profile. As 
a result, the cusped NFW profile is of no help to resolving 
the mass discrepancy between the strong lensing and X-ray 
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Table 1. Strong Lcnsing Cluster Sample 



cluster 



^cluster 



rare (Mpc) 



miens (1O 14 M ) 















fl M = 0.3 




Cl M = 0.3 


A370+ 





373 


1 


.3 


0.35 


0.41 


13.0 


13.6 











.724 


0.16 


0.19 


2.90 


2.99 


A963+ 


0. 


206 






0.0517 


0.0567 


0.25(0.21)* 


0.25(0.21)* 











.711 


0.080 


0.088 


0.60 


0.61 


A1689 


0. 


181 






0.183 


0.199 


3.6(3.0)* 


3.7(3.2)* 


A1835 


0. 


252 






0.150 


0.167 


1.98(1.54)* 


2.02(1.60)* 


A2163 


0. 


203 





.728 


0.0661 


0.0724 


0.41 


0.42 


A2218+ 


0. 


171 


1. 


.034 


0.26 


0.28 


2.70 


2.75 











.702 


0.0794 


0.0859 


0.623 


0.632 








2 


.515 


0.0848 


0.091 


0.570 


0.587 


A2219+ 


0. 


228 






0.079 


0.087 


0.517(0.415)* 


0.527(0.429)* 












0.110 


0.122 


1.60(1.28)* 


1.63(1.32)* 


A2390 


0. 


.228 





.913 


0.177 


0.196 


2.54 


2.60 


CL0024 


0. 


391 


1 


.675 


0.220 


0.257 


2.6 


2.7 


CL0500 


0. 


327 






0.15 


0.17 


1.90(1.33)* 


1.95(1.39)* 


CL2236 


0. 


552 


1. 


.116 


0.0876 


0.107 


0.30 


0.32 


CL2244 


0. 


328 


2 


.236 


0.0465 


0.0532 


0.20 


0.21 


MS0302 


0. 


423 






0.122 


0.144 


1.60(0.95)* 


1.66(1.01)* 


MS0440 


0. 


197 





.530 


0.089 


0.097 


0.89 


0.90 


MS0451 


0. 


539 






0.190 


0.230 


5.2(2.3)* 


5.4(2.4)* 


MS1008 


0. 


306 






0.26 


0.30 


6.1(4.4)* 


6.3(4.6)* 


MS1358 


0. 


329 


4 


.92 


0.121 


0.139 


0.827 


0.882 


MS1455 


0. 


257 






0.098 


0.110 


0.86(0.67)* 


0.88(0.69)* 


MS2137 


0. 


313 






0.0874 


0.0995 


0.71(0.51)* 


0.73(0.53)* 


PKS0745 


0. 


103 





.433 


0.0459 


0.0482 


0.30 


0.30 


RXJ1347 


0. 


451 





.81 


0.24 


0.28 


4.2 


4.4 



+Multipl e-arc system. 

*Arc-like image is assumed at z B = 0.8 (z s = 2). 



Table 2. X-ray properties: a standard fit to the entire X-ray surface brightness 

PLEASE PLACE TABLE 2 HERE 



measurements, and the ratios of miens /tinfw listed in Ta- 
ble 4 remain roughly the same as the values of mi e n S /m, X ray 
in Table 2. If there are any differences between the NFW 
and X-ray mass estimates, they are certainly beyond the re- 
gions accessible to current lensing and X-ray data. It could 
be misleading if one simply extrapolates our present mass 
comparison into smaller (or larger) radii. 



4 DISCUSSION AND CONCLUSIONS 

We have conducted a combined analysis of the cluster mass 
estimates from the strong lensing, the X-ray measurement 
and the NFW profile, aiming at reducing the mass discrep- 
ancy between the strong lensing and X-ray measurements 
reported in literature (see Wu & Fang 1997). In particular, 
the cusped NFW profile at the central regions of clusters was 
thought to be promising for providing a new insight into the 



problem. Basically, these three mass estimates originated 
from very different motivations, and should in principle be 
independent of each other. However, the actual application 
of the NFW profile needs the priori calibration of the two 
free parameters, p s and r s , which requires the knowledge 
about the baryonic mass distributions (galaxies and intra- 
cluster gas) in clusters, in conjunction with the spectroscopic 
measurements of the velocity dispersion/temperature of the 
clusters. So, when the X-ray data are used to fix the NFW 
profile, there is a common working hypothesis behind the 
X-ray and NFW mass estimates: i.e., the hydrostatic equi- 
librium and isothermality for the intracluster gas. Therefore, 
the gravitational lensing appears to be the unique mass es- 
timate at present independently of the dynamical state of 
clusters. 

From a detailed comparison of the projected cluster 
masses enclosed within 26 arc- like images among 21 clusters 
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Table 3. X-ray properties: cooling flow regions excluded 

PLEASE PLACE TABLE 3 HERE 



Table 4. Cluster properties: the NFW profile 



cluster a r s (Mpc) m N Fw(< r arc )* m lons /m NFW 

A370 13.62li° ^ 6 (3) 1.52±£-?°(3) 2.57;™ 5.05+f 4 20 

„ + 6!58 o V o+29.98 

U-M_ .68 d -' 8 -3.38 

A963 7.31±°-g(3) 0.31±g;g?(3) 0.15+°;^ 1.68+£;i 2 

0.29t°-S 2.10±1^ 

A1689 10.93+^^(1) 0.72±°-£(l) 1.44+° f 2.51+°,^ 
A1835 9.00+^9 (1) 0.23t°^(l) 1.44+J^ 

A2163 8.83±°i«(l) 1.03t°-°°(l) 0.3lt°-°6 1-341°^ 

+0.33 m nss+0D8 m 1/1/1+0-49 1fi7 +0.63 



A2218 9.19+^(1) 0.88tg;*(l) l-44_ . 36 



n 9 o+0.09 9 R7+0.97 

26+ 09 2 19+ 79 

A2219 11.51(2) 1.59(2) 0.341°°} l-5lt°,°6 

59+ 02 2 69+ 11 

U - Oy -0.02 Z - Dy -0.10 

A2390 8.76+°g(l) 0.54«° 3 3 (1) 1.52±°"« l-67±g-« 

CL0024 6.8ltJ° 7 *(3) 0.21«^( 3 ) 1.00^^ 2.59l°?« 

CL0500 12.91+™ 7 (3) 1.30+°;^(3) 0.73+ 9 ;^ 2.61 +22 4 °, 1 

CL2236 7.60 +2 'g(3) 0.2lt° : J 3 (3) 0.40+ 2 , , , 0.75tg : g 

CL2244 3.81±° : °j>(3) 0.22t°-' 1 1 °(3) 0.085±°-*» 2-W±lf 4 

MS0302 7.30+}; 27 -(l) 0.17+°'° 9 (1) 0.46+J'j° 3.51+^ 4 7 ° 4 

MS0440 5.62+°-^(l) 0.075+°'°}J (1) 0.31+ ,' 22 2.84±i"f| 

MS0451 5.32+°j6(l) 0.44+°;i' 4 (l) 1.01+J;^ 5.17+f;^ 

MS1008 9.03 + };^(3) 0.64 +0 ,; 2 °(3) 1-64+f J| 3.72 + |;|| 

MS1358 14.29(2) 1.48(2) ° O^t ,^} l-5llo°5 

MS1455 7.34lg;«(l) 0.13™ (1) 0.44+™° 1.96+^ 

MS2137 11.48(2) 0.18(2) °- 45 -o'u7 L59 -oi3 

PKS0745 8.93l°}°(l) 0.27 + °'° 1 1 (l) 0.23+g'^ 1-29+ ,-' 44 

RXJ1347 8.17lg;g(3) 0.18+°;° 4 (3) 2.64+ 2 ;°<; 1.59+J>* 



Nores. *In units of 1O 14 M . 

(1) From Wu & Xuc (2000). 

(2) From Ettori & Fabian (1999). 

(3) Estimatcd from a-/3 and r B -r c relations (Ettori & Fabian 1999). 



obtained from the three mass estimates, we have arrived at 
the following conclusions: 

(i) The mass discrepancy among the three methods are 
actually within a factor of 2, if X-ray measurement uncer- 
tainties are included. This is particularly true if a non-zero 
cosmological constant is invoked. 

(ii) We have confirmed the finding by Allen (1999) that 
the mass ratios of mi ens to m xray (or ttinfw) are consistent 
with unity for the massive cooling flow clusters, indicating 
that they are the more dynamically relaxed systems. 

(iii) Our X-ray mass estimates within r arc are little af- 
fected by whether or not the excess X-ray emission in the 
central cores is excluded in the fit of the observed X-ray 
surface brightness. 

(iv) The cluster masses within the arc-like images deter- 



mined from the X-ray and NFW methods show an excellent 
agreement. Any differences between these two methods must 
occur at even smaller radii (< r a rc) or at large radii. This 
arises because both methods have assumed the isothermal- 
ity and hydrostatic equilibrium for the intracluster gas. In 
other words, the employment of the cusped NFW profile for 
the dark matter distribution of clusters cannot resolve the 
reported mass discrepancy between the strong lensing and 
X-ray measurements claimed in literature. 

Taking these results as a whole, together with the previ- 
ous finding that there is good agreement between the X-ray, 
optical and weak lensing determined cluster masses on scales 
greater than the X-ray core sizes (Wu & Fang 1997; Allen 
1998; Wu et al. 1998), we feel that the reported mass dis- 
crepancy in the IC, especially NC clusters is likely to have 
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arisen from the oversimplification of the modeling of cen- Wu X.-P., Xue Y.-J., Fang L.-Z., 1999, ApJ, 524, 22 
tral mass distributions of the clusters as lenses. Indeed, if Wu X.-P., Xue Y.-J., 2000, ApJ, 529, L5 
the presence of cooling flows in the central cores of clusters 
is a reliable indicator of their dynamical relaxation (Allen 
1998), it would not be difficult to understand the reported 
mass discrepancy in the central regions of the IC and NC 
clusters where the local dynamical activities play a dominant 
role in the organization of the central matter distributions. 
Consequently, a simple spherical lens model may lead to an 
overestimate of the gravitating masses of the lensing clusters 
(e.g. Bartelmann & Steinmetz 1996; Allen 1998; Wu et al. 
1998; etc.), and meanwhile the X-ray mass estimate includ- 
ing the NFW method should also allow the intracluster gas 
to deviate from the isothermal and hydrostatic equilibrium 
distributions. 
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cluster 



Figure 1. Ratios of the strong lcnsing and X-ray determined cluster masses enclosed within 26 arc-like images among 21 clusters. The 
horizontal axis is ordered in terms of the arc list in Table 1. A flat cosmological model with (Qa = 0.7, filled squares) and without 
(Qa = 0, open squares) the cosmological constant is used. For the 11 arc-like images whose redshifts are not available, we adopt a mean 
value of z B = 2. 
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m xray ( M ©) ( entire data ) 



Figure 2. A comparison of the X-ray cluster masses enclosed within 17 arc-like images among 14 clusters (Table 3) derived from two 
different fittings of the X-ray surface brightness S x of the clusters: (1) a fi model fitted to the entire S x (horizontal axis), and (2) a (3 
model fit by excising the central region of Sx until an acceptable \ 2 is achieved (vertical axis). A fiat cosmological model of SJm = 1 is 
assumed. 



© 0000 RAS, MNRAS 000, 000-000 



10 Wu 




Figure 3. The NFW profile derived cluster masses muFw within 26 arc-like images are plotted against the corresponding X-ray masses 
(wixray)- Data are taken from Table 4. 
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Table 2: X-ray properties: a standard (3 fit to the entire X-ray surface brightness profif 



Qm — 1 Qai — 1 

cluster T (keV) ref d 

r c (Mpc) rn b x mi ens /m xray r c (Mpc) rn b x 



A370 

A963 

A1689 
Af835 
A2163 
A2218 



A2219 

A2390 

CL0024 

CL0500 

CL2236 

CL2244 

MS0302 

MS0440 

MS045f 

MS1008 

MS1358 

MSf455 

MS2f37 

PKS0745 

RXJ1347 



7 1 O + 1.05 

' • i,:J -1.05 



95+ 75 



0-48±g : g 



6.13 



+0.45 
-0.30 



9.02 
9.8 
14.6 
7.1 



0.40 
0.30 
+ 1.4 
1.4 

+0.85 
0.85 
+0.2 
-0.2 



12.4 



-0.5 
-0.5 



11 l +1 -° 

c 7+4.9 
°-'-2.1 
7 o+3.7 
' ' z -1.8 

6.2+^ 
4-6+85 

r QA+1.27 

°- ou -0.85 

in 17+ 1 - 55 



7.29±?;g 
5.45t°;» 



-0.72 

s.7±i:« 
n.37±5:g 



0.51 



+0.04 
-0.04 



65+ - 04 

U - DO -0.02 

fi^+ 04 

U - DO -0.04 

0.62+ 02 



0.65 



-0.02 
+0.08 
0.05 



40+ 07 



59+ 02 

u -^'°-0.050 

n q+°- 5 

u - a -0.3 

^+°- 18 
266+ 006 

U - ZOD -0.006 

62+ 34 

U - Dz -0.13 
-+0.03 
-0.03 
5+0.13 
-0.09 
j+0.11 
-0.07 
47+0 02 
u -^'-0.02 

64+ 04 

u -°^-0.03 



0.45] 
0.68^ 
0.63 H 



63+ 04 

u -°°-0.03 

59+0- 01 
u - oy -o.oi 

7+0.04 
-0.014 



6 
7 
1 
8 
1 
9 

10 
8 
9 
9 
9 
9 
3 



0.098 



+0.021 
-0.021 



3 





131 


4 





074 


5 





305 


5 





226 



+0.022 
-0.014 



+0.019 
-0.019 
+0.080 
-0.050 



1 47+O.O68 
u - i ^' -0.068 



14(1+0 009 

u.iiu_ 009 
066+ 039 

U.UOO_ .025 



0.41: 



0.57] 



1+0.21 
-0.24 

0.066t°;g 4 7 
0.0681°;^ 

0-075lgSS 
2 7+ - 005 

u - uz ' -0.005 

256+ 069 

U./OD_o.053 

203+ 064 

U.^UO_ .048 

042+ 013 

u - u ^ z -0.009 

62+ 009 

U.UO^_ .008 

n47+ 008 

U.IKU _ .007 

056+ 006 

U.UOD_ .006 
7+0.012 
-0.012 



n OO+4.30 
z -°°-1.56 

58+ 1 - 43 

u - oo -0.41 



0.057] 



0.13 
0.28 
1.52 
1.50 
0.22 
1.55 
0.21 
0.24 
0.33 
0.57 
1.58 
0.99 
0.57 
0.41 
0.080 
0.52 
0.36 
1.35 
1.62 
0.71 
0.51 
0.37 
0.26 
2.63 



0.05 
-0.03 
0.09 
0.06 
0.23 
0.18 
0.37 
-0.33 
+0.04 
-0.03 
+0.42 
-0.37 
+0.09 
-0.07 
+0.10 
-0.08 
0.27 
0.14 
+0.38 
0.23 
+0.25 
0.23 
+ 1.20 
0.46 
+2.07 
0.38 
+0.48 
0.21 
+0.066 
0.033 
+0.56 
0.26 
+0.12 
0.08 
+0.78 
0.49 
+ 1.16 
0.61 
+0.48 
-0.43 
+0.08 
-0.07 
0.07 
-0.09 
0.08 
0.06 
0.49 
0.30 



11.29 
3.62 
09 + 11.83 
uz -3.58 
oq+0.64 
88 -0.53 
ic+0.64 
1D -0.52 
07+O.32 
J '-0.31 
+0.37 
-0.26 
+0.31 
-0.26 
+0.53 
0.37 



5.58 H 
5 
1 
2 
2 

1.32 
1 

1.74 



0.56 



-0.44 

-0.26 



2.39 
1.59 
2 



3.35 
0.74 
2.50 
3.08 
2 

3.86 



(1.58j 



2 95+ 1 43 



+ 1.14 
-0.72 
+ 1.19 
0.72 
on+1-88 
° U -1.12 

1 61+ 027 
i - Di -0.22 

2.62+ 2 ' 31 



-1.44 

+6.67 
2.63 
+0.82 
0.40 
+ 1.76 
-1.13 
+3.01 
1.59 
+0.74 
-0.63 
+2.17 
-1.42 



53\c 
44 J 





107 


26\c 
26/ 





142 


29\c 
20/ 





083 







334 







245 



(i-28t°;|) c 
(2-25t ;^) c 



(2.34 



+4.67\ 
1.84J 



(1.83 



+ 1.79\ 
0.95/ 



-0.47 
i - ' U -0.24 

1 qn+o- 59 

i - au -0.31 
L14 +0.30 



(1.321°;?°)' 

(l-37t£g) c 



1 



0.25 
fi Q+0.21 
DU -0.25 



+0.023 
-0.023 



+0.024 
-0.015 
+0.012 
-0.012 
+0.021 
-0.021 
+0.087 
-0.054 



1fi-?+0 075 
U.10C)_o.o75 



0.155 
0.077 



0.078] 
0.088^ 



(i.7ilS:S) c 0.311 

3.771 2 ' 3 , , (2.72l 4 ^) c 0.230 

^ 7 +l.'80 



-0.010 
-0.010 
+0.045 
-0.029 

0.47;g; 24 

o.o8o±g:Sg§ 

t-0.035 
J -0.035 
h0.074 
-0.051 

29+ 005 

u.uzy_0.005 
0.084 
-0.064 
0.073 
J -0.054 

M.048tfe 

hO.010 
-0.009 
(+0.009 

i-o.oos 

o.o59j558g 
o.o68±S:Sil 



0.069^ 
0.054 H 



9 or+5.29 

z.oo_ 1 92 

74+ 1 94 
u - '^-0.53 

1 5+ 006 

u - lo -0.04 

S1+ 010 

u - ol -0.07 
-+0.25 
-0.19 
5+0.41 
-0.37 
-1+0.04 
-0.03 

1 66+ - 45 

i - DO -0.39 

23+ 010 

u - zo -0.07 



1.65] 
1.66^ 
0.24 H 



25+ 11 

u - zo -0.08 

36+ 30 

u -°°-0.15 
3+0.42 
-0.25 
-+0.27 
-0.25 
- + 1.40 
-0.54 

0.631 2 - 28 

™+°- 58 



0.63] 
1.75^ 
1.15 H 



0.091 
0.61 



+0.074 
-0.038 
+0.65 



1.62] 
1.87 H 



-0.30 

39+0- 13 

u - oy -0.09 
5+0.94 
-0.58 
7+1.33 
-0.71 

Q.8ii8:2 
o.57l8;8| 

42+ 08 
u -^ z -o.io 

n 90+0. 08 
u - zo +0.06 

3 04+ 57 

o - u ^-0.35 



a MC: massive cooling-flow; IC: intermediate cooling-flow; NC: non-cooling flow; U: Unclear. 
b In units of 1O 14 M . 

c Arc-like image is assumed at z s — 0.8 (z s = 2). 

^References: (l)Ota, Mitsuda & Fukazawa (1998); (2)Rizza et al. (1998); (3)Mohr, Mathiesen & Evrard (1999); (4)Alk 
(5)Arnaud & Evrard (1999); (6)Elbaz et al. (1995); (7)Soucail et al. (1999); (8)Schindler (1999); (9)Molikawa et al. (1999) 
(1998). 



Table 3: X-ray properties: cooling-flow regions excluded 



Q 



cluster 


T (keV) 








ref c 




r c (Mpc) 




xray 






A1689 


n no+0.40 

9-02I . 30 





T7+0 

77+ 


03 
03 


1 





97+0.05 
z 1 -0.05 


1 


9O+0.30 


2 


94 


Alooo 


Q S + 1 - 4 
9 ' 8 -1.4 


u 






n 
Z 


U 


no 


1 
1 


t-Q+0.23 
J3 -0.23 


1 

1 


9 p; 


A2163 


14 fi +0 - 85 





73t° 


02 
02 


1 





/(9+0.02 
^ z -0.02 





iq+o.03 

iy -0.02 


2 


20 


A2218 


7 -i +0.2 
' - L -0.2 





70 




2 





30 


1 


oq+o.04 

oy -0.04 


1 


95 




















16 +o.oi 

1D -0.01 

1Q+0.01 

ly -0.01 


3 
3 


81 
07 


A2219 


12 4 +0 - 5 

lz - + -0.5 





75t° 


04 
04 


3 





qqn+0.026 
°' :>y -0.026 






90+0.05 

zy -0.04 

ro + 0.10 

" JO -0.08 


1 
3 


81 
00 


A2390 


11 1+ 10 







01 
01 


4 





1 oq+0.007 
loo -0.007 


1 


,< 7+0.19 
*'-0.18 


1 


73 


MS0302 


4 6+ 08 







15 
15 


4 





nro+0.031 
uoo -0.031 





R1 +0.33 
O1 -0.26 


2 


64 


MS0440 


k c.n+1.27 

°- ou -0.85 





68« 


10 
10 


4 





nsi +0-029 
uol -0.029 





/I9+0.28 
* z -0.16 


2 


13 


MS0451 


m.ntlil 





84±° 


Of) 
00 


4 





977+0.022 
z ' ' -0.024 


1 


ro+0.49 
oo -0.36 


3 


30 


MS1008 









05 
05 


4 





9-1 q+0.024 

zlo -0.024 


1 


71 +0.86 
' i -0.51 


3 


56 


MS1358 


7.5t 4 4 ; 3 3 





83t° 


11 
10 


4 





997+0.024 
zz 1 -0.024 





ci +0.58 
Di -0.40 


1 


35 


MS1455 


K 4K+0-29 
°-^°-0.28 





66t° 


03 
03 


4 





n7n +0.007 
u(u -0.007 





ci +0.07 
dl -0.06 


1 


70 


MS2137 


^ Q7+0.38 
*- c »'-0.72 





81 




2 





09 





qo+0.03 
oo -0.06 


1 


87 


PKS0745 


8.7±i3 





65 




2 





11 





1 o+0. 03 
lo -0.02 


1 


71 



^lens / Wlxray 



r c (Mpc) 



xray 



+0.65 
-0.57 
+0.21 
-0.16 
+0.32 
-0.27 
+0.06 
-0.05 
+0.11 
-0.10 
+0.09 
-0.08 
+0.33 
-0.28 
+0.54 
-0.46 
+0.24 
-0.20 
+ 1.95 
-0.94 
+ 1.38 
-0.85 
+ 0.90 
-0.79 
+ 1.52 
-1.19 
+2.55 
-0.65 
+0.24 
-0.20 
+0.37 
-0.15 
+0.27 
-0.26 



(2.45j 
(0.97;! 



54\b 
48 ) 





9O+0.05 
zy -0.05 


1 


16\6 
127 





10 


1 







w -0.02 










32 


1 



(1.45 j 
(2.42^ 



(1.46j 
(2.57j 

(1.33^ 
(1.35j 



27 
23 


b 





375 


fO.029 
-0.029 


42 


6 








39, 

















202 


f0.008 
-0.008 









068 


f0.036 

-0.036 









089 


f0.032 
-0.032 


44 
35 







336 


f0.027 
-0.029 


10 
86 







242 


fO.027 
-0.027 









260 


fO.027 
-0.027 


19 
16 


6 





n78+0- 008 
u,o -0.008 


27 
11 


6 






10 
12 





0.21: 
0.31: 



0.71] 
0.45j 



0.18] 



°In units of 10 14 M Q . 

6 Arc-like image is assumed at z s = 0.8 (z s = 2). 
c Rcfcrences: (l)Vikhlinin, Forman & Jones (1999); 



(2)Ettori & Fabian (1999); (3)Rizza et al. (1998); (4)Lewis et al. 



